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model treated in an inhomogeneous unrestricted Hartree-Fock approximation, we find that corre- 
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which causes the resistivity upturn. A pseudogap in the density of states is shown to enhance both 
the disorder-induced magnetic state and the resistivity upturns. 
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I. INTRODUCTION 

Exotic transport properties in metallic hole-doped 
cuprates reflect their strongly correlated nature over a 
large part of their phase diagram. In the vicinity of opti- 
mal doping, the in-plane resistivity is found to be linear 
in temperature T, with deviations from linear-T power 
laws evolving on the underdoped and overdoped sidesii*^ 

\ As the transition temperature Tc is suppressed down to 
zero by a magnetic field, a resistivity that diverges log- 
arithmically at low temperatures (log-T) is observed in 
La2-a;Sr2;Cu04 (LSCO)^i^ across a wide range of dop- 
ing. Resistivity "upturns", increasing p as temperature 
T is decreased below a temperature T,„i„ have been ob- 

' served as well in Ba2Sr2_a; La^rCuOg+A- (BSLCO)^i^ and 
sufficiently disordered and underdoped YBa2Cu307_5 
(YBCO) samples.-i'^'aiio>iiii2a3ii4.i5i6 guch upturns are 

frequently associated with a metal-insulator transition as 
the system approaches its antiferrornagnetic(AF) parent 
compound (for a review, see Ref. llTI ). However, one 
should keep in mind that the explanation of these resis- 
tivity upturns must include not only the intrinsic elec- 
\ tronic correlations present in the system, but also their 
interplay with the external perturbations introduced to 
suppress superconductivity. 

\ Besides the suppression of it is known from in- 
elastic neutron scattering (INS), nuclear magnetic res- 

' onance (NMR), and muon spin rotation (//SR) experi- 
ments that introducing disorder and magnetic field can 
induce local magnetic order, reflecting the coexistence 
of strong AF correlations with superconductivity. For 
instance, a strong signal centering at incommensurate 
positions near (tt, tt) has been observed in INS experi- 
ments in the presence of a magnetic fiel d^^i^^'^°i^^ , in- 
dicating the formation of AF order around vortices in 
LSCO; a smaller but significant signal is also present in 
zero field. Other neutron scattering measurements have 
observed evidence of ordered static magnetism in intrin- 
sically disordered cuprates, and shown that systematic 



addition of disorder enhances this effect i^^'^'^i^ '^^ NMR 
measurements have furthermore shown evidence that lo- 
cal magnetic moments are induced around atomic scale 
defects such as Zn substitutions of planar Cu, or defects 
produced by electron irradiatiou i^'^'^^i^^'^^i^^ The suscep- 
tibility of these induced moments shows a Curie- Weiss 
behavior even though the impurity itself is nonmagnetic, 
indicating their origin in the strong magnetic correla- 
tions present in the pure system. Finally, /iSR experi- 
ments have shown that the Cu spins freeze in the un- 
derdoped superconducting state, and eventually develop 
short range order at very low temperatures in intrinsi- 
cally disordered cuprates and even in the much cleaner 
system YBCO if it is highly underdoped,^i^i^i23i34i35 
The relationship between ordinary disorder and local 
magnetism in these and other experiments, has been re- 
viewed in Ref. together with a description of recent 
theoretical work. Since these phenomena are well estab- 
lished, a theory which seeks to account for the transport 
anomalies should therefore also be capable of explaining 
the formation of these local moments, as well as their 
ordering behavior at different dopings and temperatures. 

The logarithmic temperature dependence of the resis- 
tivity upturns in a magnetic field has remained a mystery. 
It is tempting to associate these logs with the quantum 
corrections to the conductivity found in weak localiza- 
tion theorjj^i^ii^. Indeed, in electron-doped cuprates'^-, 
where interaction effects are thought to be weaker and 
disorder effects stronger, as well as in overdoped cuprate 
samplesi^, good fits of the magnetoresistance data to 
weak localization theory have been obtained. By con- 
trast, elastic free paths in hole-doped samples are much 
larger than the Fermi wavelength scale required for weak 
localization effects; furthermore the magnetoresistance 
has the wrong field dependence and typically (but not 
always^) the wrong sign. A log-T behavior of the resis- 
tivity is also found in the theory of granular systems;^ 
but evidence for granularity in the conventional sense is 
weak or absent in the cuprate samples where the up- 
turns have been observed. Finally, it has been argued by 
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Alloul and others that the body of experimental results 
on underdoped cuprates, specifically Zn-substituted and 
irradiation damaged YBCO samples, is consistent with 
a one-impurity Kondo picture, with conventional resis- 
tivity minimum. However there are several inconsisten- 
cies associated with this approach, reviewed in Ref. 
We adopt here the alternate point of view that the up- 
turns observed in the underdoped, hole-doped cuprates 
are manifestations of disorder in a Fermi liquid in the 
presence of strong antiferromagnetic correlations. 

A theory that can cover the anomalies of transport 
properties in cuprates over a wide range of doping does 
not currently exist. Recently, an attempt was made to 
treat disorder and interactions in a model tailored to the 
cuprates by Kontani et alM-t^^^ Within the fluctuation- 
exchange (FLEX) approach and certain approximations 
regarding the impurity scattering processes and self- 
energy, these authors had considerable success in repro- 
ducing resistivity upturns observed in some cuprates in 
zero magnetic field. However, as a perturbative approxi- 
mation it perforce neglects certain self-energy and vertex 
correction diagrams; in addition, the physical content of 
the approximations made is not always clear. 

Here we focus on the optimally- and slightly under- 
doped cuprates, in the spirit of Kontani et ali^, and 
assume the Fermi liquid picture properly describes the 
electronic excitations in the normal state. We examine 
the following simple hypothesis that connects the trans- 
port anomalies with the impurity induced magnetization: 
the resistivity upturns are due to the extra scattering as- 
sociated with the correlation-induced magnetic droplets 
which carry local moments. Within a 2D single band 
Hubbard model where interactions are treated in mean 
field but disorder is treated exactly, we show that the re- 
sistivity increases coincide with the conditions which en- 
hance impurity induced magnetic moments. The present 
study focusses on the doping regime where static mo- 
ments, even in most strongly correlated LSCO, are para- 
magnetic centers induced by the applied field. Other re- 
cent studies relevant to this phase have examined the 
more disordered, or more correlated state where such 
magnetic droplets are spontaneously formed around de- 
fects in zero field, and shown that they can indeed af- 
fect macroscopic observables such as NMR, thermal con- 
ductivity and superfluid density k{T) M^^^^^^^^L^^ The 
physical picture of the ground state, that of an inhomo- 
geneous mixture of AF droplets carrying net moments 
near the defect, is quite similar in our case. By working 
in the regime where moments are smaller and the effect 
of the field is larger, however, we hope to explain some of 
the observed puzzling aspects of the magnetoresistance. 
Since we consider relatively weak correlations, we explic- 
itly confine ourselves to the doping regions in each system 
under consideration where the resistivity upturns first set 
in. This means that we, within a RPA treatment of the 
correlations, do not expect to be able to describe the 
true MIT or log-T behavior, but rather the leading per- 
turbative corrections to the high-T behavior of the re- 



sistivity. The conditions in which positive correlations 
between impurity-induced magnetization and transport 
anomalies can be found are examined, which confirm our 
hypothesis that the enhancement of the scattering rate 
is due to an enlarged cross sections associated with these 
induced moments. We first examine the case of optimal 
doping, and then discuss the effect of including a pseudo- 
gap in the density of states, which will allow us to extend 
the model to lower dopings. 



II. MODEL HAMILTONIAN 

Since resistivity upturns are revealed after Tc is sup- 
pressed down to zero, the pairing correlation is ignored in 
describing the normal state properties as a first approx- 
imation. We therefore start with the two dimensional 
Hubbard model to describe the Cu02 plane 

H ^^-tijcl^Cja + ^{eia-f^)nia + ^Uni'ihii , (1) 



where Ci^ is the electron operator at site i with spin cr, 

flier = c|^Cicrj tij = t, t' is the hopping amplitude between 
nearest-neighbor {t) and next-nearest-neighbor (<') sites, 
and U is the onsite Coulomb repulsion. The external per- 
turbation due to impurities and magnetic field is included 
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(2) 



where Vr is the scattering potential produced by defects 
such as Zn substitution or electronic irradiation. The 
Zeeman term takes into account the spin-dependent en- 
ergy shifts caused by the magnetic field with a = +/ — for 
spin up/down respectively. We will denote ^g^sB = B 
in the figures presented below. A Hartree-Fock mean field 
decomposition is then adopted to the above Hamiltonian 



m,. 



(3) 



and gives rise to 



ij(7 ia 

U „ ' n,a ■ (4) 



By applying this mean-field Ansatz we certainly cannot 
study Mott correlations in the regime U ^ t. We can, 
however, discuss qualitatively the physics arising from 
the tendency to form magnetic moments near impuri- 
ties, as is appropriate near optimal doping. To study the 
response of the system in a static electric field, we cal- 
culate conductivity via linear response theory, where we 
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use the current in the x-direction, 

a 

+ ^(cl+2;±yo-^i'T ~ ^L^i+2;±i/o-) , (5) 
a 

and express the site-dependent current-current correla- 
tion function in terms of eigenstates and eigenenergies 



^,,{t) = -iQ{t){[Mt),j,m) , 

■nij{uj) = {n\J.i\m) (to| Jj|n) 



UJ + En - Em + ir] ' 



= y - hm 



Im (TTy (tj)) 



= 7r^^(n|J,|m) (m|Jj|n)F(£;„,£;™) . (6) 

j n.m 

The global conductivity a is realized by averaging ai 
over the whole sample with a proper normalization. 
The function F{En,Em) is symmetric under exchange 
of En <-> Em ^ The resistivity p is then given by the in- 
verse of tJ, and is plotted in units of 2D resistivity fi/e^. 
One can also convert it into a 3D resistivity for materi- 
als such as YBCO, in which one assumes two conducting 
planes per unit cell and gives p as 3D resistivity in units 
of 241/iri-cm. Note that this procedure gives us only 
the resistivity part due to impurity scattering; since the 
Hamiltonian ^ is decoupled at the Hartree level, the in- 
elastic processes which lead to, e.g. the linear resistivity 
at optimal doping are not treated. We calculate therefore 
only the low-T part due to elastic scattering. 

The proper choice of system size in simulating Eq. ^ 
is determined by the following criteria. Firstly, in the 
absence of impurities, homogeneous resistivity po should 
be proportional to the artificial broadening 77. Secondly, 
the resistivity in the case with impurities should be pro- 
portional to the impurity concentration. We found that 
a 40 X 40 lattice is able to achieve the above two crite- 
ria down to temperature as low as Ts = 0.02, which is 
roughly equal to the average energy level spacing, and 
will be the system size used in the following. Each data 
point is then averaged over 10 different impurity configu- 
rations, which we found to be sufficient to ensure the ran- 
domness of the impurity distribution. We chose t' — —0.2 
and the energy unit to be t = lOOmeV, which gives tem- 
perature scale T = 0.0 It ~ lOK and the magnetic field 
scale B ~ 0.004 ^ 7T, in the same scale as a recent study 
of NMR lineshapes in the superconducting statei^ 

Before the resistivity under the influence of induced 
magnetization is studied, we first compare the present 
study in the normal state with the data in the d-wave 
superconducting (dSC) statei^ Such a comparison re- 
veals the importance of finite DOS in the normal state, 
as well as the bound state formation in the dSC state. 
The effect of nonmagnetic impurities in the dSC state is 
studied within the framework of d-wave ECS theory plus 



magnetic correlations, equivalent to the Hubbard model 
in Eq. ([1]) with additional pairing correlations between 

nearest-neighbor sites Hpatr = Yl,iS ^Si'^\-\'^\+5i + H.c, 
where the gap is to be determined self-consistently A^i = 
V {ci^Ci+si + Ci+5|Ci|)/2, with V = 1. The real-space 
magnetization pattern induced by a single nonmagnetic 
impurity is shown in Fig. [IJa) and Fig. [IJb), in which 
we found three major differences. (1) In the presence 
of a magnetic field, the normal state has homogeneous 
magnetization significantly larger than that of the dSC 
state. This is obviously due to the opening of the gap in 
the dSC state which reduces the DOS at the Fermi level, 
and hence exhibits a smaller susceptibility than the nor- 
mal state. (2) The magnetization on the nearest-neighbor 
sites of the impurity is drastically enhanced in the dSC 
state, which we found to be consistent with the bound 
state formation due to the c?-wave symmetryJi (3) The 
dSC state has a shorter correlation length, resulting from 
the enhancement of nearest-neighbor site magnetization 
in comparison with the relatively smaller magnetization 
on the second and third nearest sites away from the im- 
purity. To give a quantitative description of these fea- 
tures, we introduce the total magnetization Sz and the 
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FIG. 1: (Color online) (a,b) Real space magnetization pat- 
tern induced by a single nonmagnetic impurity for the nor- 
mal state (a), and the dSC state (b), both at f/ = 1.75 and 
B = 0.01. One sees that the dSC state has more pronounced 
nearest-neighbor site magnetization due to bound state for- 
mation, and has smaller homogeneous magnetization due to 
the opening of a gap at the Fermi surface. These effects are 
shown more clearly in the total magnetization Sz (c) and the 
magnetic contrast A (d) versus external field. 
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magnetic contrast A 



A= -^^\m,-mo\ 



homogeneous magnetization in the absence of impurities 
but in the presence of a magnetic field. The meaning of 
A is to estimate the fluctuation of site-dependent mag- 
netization away from its homogeneous value mo, hence 
an indication of locally induced staggered moment. Since 
interference between impurities is always present and the 
local environment is different around each impurity, the 
deviation from mo of the whole system needs to be con- 
sidered, and therefore we sum over i for A in Eq. ([7]). The 
behavior of Sz and A versus the applied field is shown in 
Fig- HJc) and Fig. [Hd), where one sees that Sz in the 
normal state is one order of magnitude larger than in the 
dSC state, which is attributed to the overall larger ho- 
mogeneous magnetization in the normal state. However, 
in the A versus field plot, we see that after the homoge- 
neous magnetization is subtracted, as in the definition of 
A, the dSC state has a larger value due to the enhanced 
magnetization attributed to the bound state formation. 
Such a comparison indicates that DOS at the Fermi level 
is crucial to the formation of impurity induced moments, 
which in turn motivates us to propose a phenomenologi- 
cal model that emphasizes the effect of reducing DOS in 
the underdoped region, as will be discussed in Sec IV. 



III. RESISTIVITY UPTURNS AT OPTIMAL 
DOPING 

Motivated by the NMR experiments,iii2£i21i2ai2a we 

study the magnetic response in the paramagnetic region 
close to the magnetic phase boundary. For convenience 
and direct comparison to experiments where unitary scat- 
terers are created by Zn substitution or irradiation de- 
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FIG. 2: (a) Magnetic contrast and (b) resistivity versus U at 
optimal doping with T = 0.03 and 2% impurities. 
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FIG. 3: (a) Magnetic contrast A and (b) change of resistivity 
Ap/po versus T at optimal doping with U = 1.74 and 2% 
impurities. 



contrast A versus U. For the band structure used in 
this paper, the critical Coulomb repulsion is found to 
be Uc ^ 1.75, above which a spontaneous magnetization 
is observed for zero field. This value is found to depend 
on system size and impurity content, but the value is 
roughly close to Uc ^ 1.75. As seen in Fig. [2Jb),the resis- 
tivity increases with U, and coincides with the behavior 
of A in the region both below and above its critical value. 
This positive correlation between A and p serves as the 
first evidence that we can attribute the increase of re- 
sistivity to the extra scattering induced by the magnetic 
moments. In the following discussion we choose U = 1.74 
such that it is close to but slightly below the critical Uc, 
and the system exhibits paramagnetic response to an ex- 
ternal field. 

We note that in the region where the system cross 
the magnetic phase boundary, for instance at large U 
or low temperatures, numerics found that there are sev- 
eral stable states with comparable energies competing 
with each other. Taking different initial conditions or a 
different route for the convergence can result in a dif- 
ferent apparent ground state configuration; for instance, 
we found a charge density wave (CDW) ground state 
with periodicity (tt, tt) that can exist in large U and zero 
field, consistent with the spin or charge modulated state 
found in other studies with a sufficiently large Coulomb 
repulsion4^i^iSi^i^ However, considering the strong 
experimental evidence of magnetic ordering, as well as 
the resulting resistivity in comparison with the trans- 
port measurement, only the paramagnetic induced mo- 
ment state can give a proper description of both induced 
magnetization and transport anomalies in the optimal to 
lightly underdoped systems, and hence will be the stable 
configuration focused on in this report. 

Due to the limited system size, we are unable to ex- 
plore the extremely low T regime, which prevents us from 
comparing the present theory with the experimentally 
observed Log-T divergence. However, numerics down to 
as low as T = 0.026 ^ 26K shows significant resistivity 
upturns in comparison with the zero field case. Figure [3] 
shows both magnetic contrast and change of resistivity 
Ap/po versus temperature T, where po is the resistivity 
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at the uncorrelated zero field case(C/ = 0,B = 0), and 
one sees again the positive correlation between these two 
quantities. The lowest temperature explored is slightly 
lower than the critical temperature Tspon ~ 0.025 below 
which a spontaneous magnetization is observed in the 
zero field. The magnitude of the upturn at T = 0.026 
in comparison with high temperature resistivity is of the 
order of 5%, roughly consistent with the value obtained 
in slightly underdoped YBCO after the linear-T contri- 
bution has been subtracted^, 

The magnetoresistance in the presence of induced mag- 
netization is shown in Fig. |4l where we again see a pos- 
itive correlation between A and Ap/po with increasing 
magnetic field B. At the temperatures where the resis- 
tivity upturns set in, we found that both A and Ap/po 
first increase with the field, and eventually saturate and 
slightly decrease in the high field region. One can unam- 
biguously define a field scale B^at above which A and 
Ap/po saturate, and we found that B^at decreases as 
temperature is lowered. Such a increase-saturation be- 
havior is consistent with the magnetoresistance observed 
in YBCO,^^'^- although Bgat observed therein is slightly 
higher, possibly due to the higher field required to elim- 
inate the superconductivity before normal state proper- 
ties can be observed. Since Bgat decreases as lowering 
temperatures, the region where the magnetic contrast A 
is linear with respect to the external field also decreases 
accordingly, which indicates that as the magnetization 
starts to grow at low temperatures, the interference be- 
tween the magnetic islands induced around each impurity 
is also enhanced, causing A to deviate from a linear re- 
sponse. 

The last issue we need to address is the behavior of A 
and Ap/po as changing impurity concentration riimp, in 
comparison with the available experimental data which 
shows that the resistivity upturns monotonically increase 



with Tiimp up to rii. 



3%. Fig. [5] shows the numer- 



ical result under the influence of changing riimp, where 
one again sees the consistency between the behavior of 
A and Ap/po. However, instead of increasing monoton- 
ically with increasing Uimp, we found that both A and 0.03 
Ap/po increase up to a critical concentration n: 



imp 



1%, 

and then decrease as more impurities are introduced on 
the plane. Such a result indicates that the impurity in- 
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FIG. 4; (a) Magnetic contrast and (b) change of resistivity 
versus B at optimal doping with U = 1.74 and 2% impurities. 



duced magnetization is proportional to riimp only up to 
a certain extend, beyond which the interference takes 
place and eventually destroys the magnetization and the 
associated magnetic scattering. To further demonstrate 
that the interference effect is more destructive than con- 
structive to the induced magnetization, we study the 2- 
impurity case in the present model, and plot A against 
the separation between the two impurities, as shown in 
Fig. [DJb). We first found that there exists a strong en- 
hancement of magnetization if both impurities are on the 
same sublattice, consistent with previous studies in the 
dSC state4^i^ Secondly, A indeed decreases as the two 
impurities get closer, which is the case when ntmp is in- 
creased, indicating the destructive nature of the interfer- 
ence effect, and hence the decreasing of magnetization 
at sufficiently large impurity content. Our result there- 
fore predicts that if the extremely disordered samples 
{nimp > 3%) can be studied experimentally, a critical 
concentration can occur beyond which the resistivity up- 
turn drops as increasing impurity content, assuming that 
weak localization has not yet taken place. The critical 
concentration Uimp ~ 1% shown in the present study is 
apparently smaller than the experimental value, which 
may be due to a smaller linear response region in the 
present model in comparison with the real cuprates, pre- 
sumably an artifact of such a weak coupling mean field 
approach. In addition, the critical disorder concentration 
TT-imp will depend on the details of the disorder modeling, 
for instance the nature of the disorder, or the extent of 
the impurity potential, which is outside of the scope of 
our study. 



IV. EFFECT OF PSEUDOGAP IN DOS ON 
RESISTIVITY 



From a weak coupling perspective, we expect enhance- 
ment of resistivity upturns as the system is underdoped. 
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FIG. 5: (a) Magnetic contrast A and change of resistivity 
Ap/po versus riimp at optimal doping with U = 1.74, B = 
0.001, T = 0.03, and (b)A induced by the 2-impurity model 
plotted against the seperation between the two impurities ri2 , 
collecting all relative positions up to thirdteenth shell. Values 
of ri2 that correspond to average distance of impurities at 
riimp ~ 3%, 1%, and 0.5% are indicated. 
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based on the following two features: Firstly, correlations 
are more prominent as one goes toward half-filling, result- 
ing in an increase of the effective U entering our model. 
Although the Hartree-Fock type mean field theory can 
not capture the Mott transition induced by correlations 
nor the pseudogap phenomenon, the drastic increase of 
resistivity near the critical value of U suggests that cor- 
relations indeed affect resistivity as one approaches the 
strong coupling region. The large U region in Fig. [2] 
demonstrates that correlation strength J7, as well as the 
induced magnetic moment, are indeed essential ingredi- 
ents to determine the magnitude of the upturn. 



DOS is proposed for the homogeneous pseudogap state 



El, 



N{uj) = 



sign(Cfe)\/^fc 



dp 



(8) 



where ^fc ~ —2t{cos{kx) + cos{ky)) — At' cos{kx) cos(fcy) — 
/i/ is the normal metallic dispersion, with a constant 
"pseudogap" = 0.2. We then Fourier transform Ek 
back to real space and find an effective long range hop- 
ping model that gives the energies E^ ■ The hopping am- 
plitude tij of this extended hopping model is therefore 
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FIG. 6: Comparison of DOS: the normal state dispersion ^fc 
gives NiiuS) (shifted), the proposed phenomenological model 
for the pseudogap state gives N2{uj) (shifted), and the 
actual A^3(a') (original scale) given by the effective hopping 
model after Fourier transform of Ek in a 40 x 40 system, with 
Hf = 0.02 and A = 0.2. 



Secondly, the opening of the pseudogap in the quasi- 
particle spectrum is known to favor bound state forma- 
tion, which in turn promotes the impurity induced mag- 
netic moment.~ This is similar to the dSC state where 
the pole of impurity T-matrix falls within the gap, pro- 
ducing a bound state localized around the impurity. We 
expect that the reduction of the DOS in the pseudogap 
state also produces poles of T-matrix near Fermi energy, 
although the exact form of Green's function and Dyson's 
equation remains unknown. Resistivity upturns are then 
affected by the pseudogap formation, based on the naive 
argument that impurity induced moments result in the 
upturn. To get a crude idea of the effect of reducing the 
DOS, we introduce a pseudogap in an ad hoc way without 
going through the T-matrix formalism, since no micro- 
scopic model of the pseudogap state is generally agreed 
upon at present. The following form of dispersion and 



4^ 



Ek{cos[kx-{xi-Xj)]+cos[ky-{yi-yj)]} . (9) 



We calculate the hopping range up to \xi — Xj\ = 
lUi — Ujl = 20 on a 40 X 40 lattice. Numerics show a 
roughly 40% reduction of DOS at the chemical potential, 
as shown in Fig. [S] The calculation of resistivity then 
follows Eq. ([5]) and while the contribution from all 
hopping terms ts = Uj and their corresponding distance 
d = fi — fj all need to be considered. 
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(Extended) and without (Normal) reduction of DOS by apply- 
ing extended hopping Eq. (|9]), both at optimal doping, with 
U = 1.75, B — 0.001 and 2% impurities, (c) and (d): same 
quantities vs. B for Normal and Extended hopping models. 

Figure [7] (a,b) shows the magnetization and resistiv- 
ity comparing extended hopping model with the normal 
state Hubbard model Eq. ([T]) which contains only nearest 
and next-nearest neighbor hopping. Wc fix both models 
at optimal doping 6 = 0.15 and examine solely the ef- 
fect of reducing DOS. Among the magnetic field region 
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explored < < 0.01, the magnetic contrast A is 

found to be enhanced in the extended hopping model, 
confirming our hypothesis of reducing DOS promotes 
bound state formation, which also gives slightly larger 
resistance between temperature range 0.02 < T < 0.045. 
The magnetization and resistivity versus field is shown in 
Fig. [7Jc,d), where one sees larger magnetization compar- 
ing to the normal state model, with a smaller linear re- 
sponse region and the saturation at high field is again re- 
vealed. Resistivity upturns are enhanced overall in both 
low and high field region, and is consistent with the be- 
havior of A. The hypothesis of reducing DOS promotes 
induced moments, and in turn enhances the resistivity 
upturns, is then well proved. 

V. CONCLUSIONS 

In summary, we employed a Hartree-Fock decompo- 
sition of the Hubbard model to study transport prop- 
erties under the infiuence of disorder induced magneti- 
zation, which is a consequence of the interplay between 
strong correlations and inhomogeneity. The numerical 
results suggest that, at low enough temperatures and 
strong enough correlations, impurity induced magneti- 
zation is drastically enhanced. Within this regime, both 
induced magnetization and resistivity are increased as 
(1) the temperature is lowered, (2) the magnetic corre- 
lations are enhanced, (3) the magnetic field is increased, 
and (4) more impurities arc introduced, consistent with 
the conditions in which the enhancement of resistivity is 
observed experimentally. We predict, in addition, that 
the addition of further disorder can sometimes lead to a 
nonmonotonic field dependence as the magnetic poten- 
tial landscape becomes smooth; this property has not 
yet been observed to our knowledge. Extremely heav- 
ily disordered or strongly correlated samples will lie in a 
different regime, which we have not yet treated, where 
disorder will create a spontaneous, short-range ordered 



magnetic state even in zero field*^ ; in this case we antic- 
ipate that the magnetoresistance will quite small. 

The positive correlation between induced magnetiza- 
tion and resistivity confirms our hypothesis that the en- 
larged cross section due to these local magnetic moments 
gives extra scattering and hence the resistivity upturns, 
and indicates that the hole-doped cuprates lie within this 
regime over a wide range of (under) doping, in which 
strong correlations can cause anomalies in the thermo- 
dynamic observables. A phenomenological model that 
produces reduction of DOS near the Fermi level in an 
ad hoc way further suggests that, as the system is un- 
derdoped, besides the enhancement of correlations that 
can increase the resistivity, the anomalous energy spec- 
trum in the underdoped region can promote the impurity 
bound state and hence the magnetization, which in turn 
boosts the magnetic scattering and the resistivity up- 
turns. The proposed mean field theory plus real space 
diagonalization scheme is therefore a powerful tool to 
capture the complex effect on the transport properties 
due to strong correlations, inhomogeneity, and the spec- 
tral anomalies in the low temperature region where the 
transport is dominated by disorder. Further applications 
of the present theory, as well as the influence of impu- 
rity induced magnetization on other thermodynamic ob- 
servables in the metallic cuprates, will be addressed in a 
future study. 
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